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METHOD FOR FINDING OPTIMAL THRESHOLD 
FOR IMAGE SEGMENTATION 

Field of the Invention 

5 

The present invention relates to a method for finding a 
threshold for image segmentation; and, more particularly, to a 
method for finding the optimal threshold for image 
segmentation in image recognition. 

10 

Description of Related Arts 

Generally, a process of finding the optimal threshold 
for the image segmentation is fundamental and important 
15 process in the image recognition . The process of recognition 
is necessary to distinguish an object from a background of an 
image . 

The optimal threshold can be found easily based on a 
bimodal type histogram distribution graph and in above case, 

20 it is located at a lowest point of histogram distribution 
curve. There are many methods introduced for finding the 
optimal threshold. 

A first method is stochastic method to find the optimal 
threshold. That is, a histogram distribution of an image is 

25 assumed as the bimodal type and a gray level having the 
minimum sum of variance is select as the optimal threshold. A 
second method finds the optimal, threshold based on the Shannon 
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entropy. A gray level having the minimum entropy is selected 
as the optimal threshold. There is also another method using 
fuzziness during computing entropy of image. This method 
selects a gray level having minimum fuzziness as the optimal 
5 threshold. 

Fig. 1 is a flowchart explaining a conventional method 
for finding the optimal threshold. Referring to Fig. 1, at 
step S101, a histogram distribution of an image is computed. 
Entropies of all gray levels distributed in the histogram are 

10 computed at step S102. All computed entropies are compared 
one another and a gray level having the lowest entropy is 
selected at step 103. However, a process time of the 
conventional method is increased in proportion to a resolution 
of image and the number of entropies to be computed. 

15 Furthermore, in case of an image having multi thresholds, a 
* condition of segmentation completion may not be computed by 
the conventional method. 

Summary of the Invention 

20 

It is, therefore, an object of the present invention to 
provide a method for effectively finding an optimal threshold 
for image segmentation of an image having multi thresholds by 
analyzing entropy characteristic of image based on a fixed 
25 point iteration and fuzzy entropy. 

In accordance with an aspect of the present invention, 
there is provided a method for finding a threshold value in 



image segmentation, the method including the steps of: a) 
gaining histogram distribution of an image; b) computing 
entropy values corresponding to gray levels in the histogram; 
and c) gaining a minimum entropy value corresponding to the 
5 gray level as the threshold value by using a fixed point 
iteration FPI based on the computed entropy values. 

Brief Description of the Drawing (s) 

10 The above and other objects and features of the present 

invention will become apparent from the following description 
of the preferred embodiments given in conjunction with the 
accompanying drawings, in which: 

Fig. 1 is a flowchart explaining a conventional method 

15 for finding the optimal threshold; 

Fig. 2 is a flowchart for explaining a method for 
finding the optimal threshold for image segmentation in 
accordance with a preferred embodiment of the present 
invention; 

20 Fig. 3 is a graph showing a fuzziness according to a 

constant k in accordance with a preferred embodiment of the 
present invention; 

Fig. 4 is a graph showing gray level distribution curve 
for obtaining a gray level of minimum entropy in accordance 
25 with a preferred embodiment of the present invention; 

Fig. 5 is a flowchart for explaining the step S203 in 
Fig. 2 for gaining a gray level corresponding to the minimum 
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entropy by using FPI in accordance with a preferred embodiment 
of the present invention; and 

Fig. 6 is a flowchart for explaining step 501 of Fig. 5 
in detail. 

5 

Detailed Description of the Invention 

Other objects and aspects of the invention will become 
apparent from the following description of the embodiments 
10 with reference to the accompanying drawings, which is set 
forth hereinafter . 

Fig. 2 is a flowchart for explaining a method for 
finding the optimal threshold for image segmentation in 
accordance with a preferred embodiment of the present 
15 invention. 

Referring to Fig. 2, a histogram distribution of an 
image is obtained at step 201. Entropies of gray levels are 
computed at step S202. After computing entropy value at step 
S202, a gray level of minimum entropy is gained by using the 
20 fixed point iteration (FPI) based on the computed entropy 
value at step S203. 

In the step S202, the entropy of the gray level is 
computed by measuring fuzzy entropy of corresponding gray 
level. Hereinafter, the computation of fuzzy entropy is 
25 explained in detail. 

If there is an M x N size of image I having L gray 
levels, a gray level of pixel (x,y) is defined as I(x,y) and 
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JUjUij) represents fuzziness of gray scale of pixel (x,y). 
Therefore, the image / can be expressed as following equation. 

/ = {(/,,/!,(/,))> Ec *' 1 

5 , wherein 0 < //,(/,,,) S 1;/ = 0,1,..., M -1; j = 0,1,..., AT- 1 

If a gray level g has a frequency of generation h(g) in 
entire image / then an average gray level // 0 of a background 
can be expressed as following equation 2 and an average gray 
10 level Hi of an object can be expressed as following equation 3. 

M>=i>Ks)/l>(s); if Kx,y)ZT Eq. 2 

ft=I^)/&(g); if i(x,^)<r Eq. 3 

g=r+i r+i=o 

15 

The average gray levels n l and fi Q can be considered as 
target values of the threshold value T. That is, the 
fuzziness can be expressed as a difference between gray level 
of a pixel (x,y) and a gray level of a region including the 
20 corresponding gray level. Therefore, the difference of gray 
levels is smaller, as larger as the fuzziness is. Gray levels 
of all pixels in an image for a threshold T must have certain 
fuzziness either to an object or background. The fuzziness of 
a pixel can be expressed as following equation. 
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! if i(x,y)<T 

\+\i(x,y)-Mo\lk e ^ 4 

if i(x, > ar 

1+|i(x,jO-/jJ/* 

In a meantime, when a gray level of a certain pixel is 
5 included in a specific region, the fuzziness must to be more 
than 0.5 to a corresponding set. Therefore, a constant K is 
defined as G max -G min and the fuzziness is selected between 0.5 
to 1. Fig. 3 shows fuzziness according to the constant K. 

Also, fuzziness can be generally considered as a 
10 fuzziness of fuzzy set and there are various entropy-equations 
disclosed for computing fuzzy entropy. If the entropy 
equation of one independent variable is expanded to 2 
dimensional image region then entropies of the image I can be 
expressed as following table 1. 

15 

table. 1 



Entropy 
type 


Entropy Measures 


Shannon 


MNin(2) V , 
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Yager 


s<« (<(*■ rt)) = m, ('(*. ^)) - y)) 

wnprp <y i<s a fii77ifipr factor 


Pal&Pal ! 


= Y Y (i(x, y))), 

S(u, (i(x,y))) = u. (i(x, y))»exp[ u, (i(x, y))l 

+ { [1 + 0' (*> ^))3-exp[l + //, (/(x, y))] } 


Bhandari 


£(./) = 1 y y (/(x, y))\ 

MVln(2)(l-a) V v 
S( Ml (/(*, y))) = log[//, (/(x, y ))" + (1 - Ml (,(x, y)) a ] 


Standard 
tuzzy 

Complement 


M-LV X y 

S(m, (/(x, y))) = 1 - [2 M , (i(x, y)) - 1] 


Kaufmann 


^ J > = t^EZ 5 ^/ (*(*.*))), 

S(M, ('(*, y))) = min{// / (i(x, y)), 1 - //, (/(x, y))} 


Quadratic 
Kaufmann 


f T 1/2 
S(v, (.Kx, y))) = mini//, (/(x, y)), 1 - //, (/(x, y))} 2 



Specially, the absolute value of entropy is increased in 
a range of [0, 0.5] and the value of entropy is decreased in a 
region of [0.5, 1]. In a meantime, if fuzziness of all gray 
level included in the image are about 0.5 then entropy E(I) 
has 1 as the maximum value. 

After obtaining entropy values at step S202, a gray 
level of minimum entropy is obtained by using a fixed point 
iteration (FPI) method at step S203. 



Fig. 4 is a graph showing gray level distribution curve 
for obtaining a gray level of minimum entropy in accordance 
with a preferred embodiment of the present invention. 

Referring to Fig. 4, obtaining a gray level of minimum 
5 entropy is explained hereinafter. 

Fig. 5 is a flowchart for explaining the step S203 in 
Fig. 2 for gaining a gray level corresponding to the minimum 
entropy by using FPI in accordance with a preferred embodiment 
of the present invention . 
10 At step 501, possible optimal threshold values P± for 

obtaining a gray level of minimum entropy are obtained based 
on the graph of Fig. 4. 

After obtaining the Pi, gray levels are sequentially 
obtained from left to right Pi. 
15 At step 503, an optimal threshold of gray level having 

minimum entropy is obtained by comparing entropy values of 
gray levels of Pi. 

Fig. 6 is a flowchart for explaining step 501 of Fig. 5 
in detail. 

20 Referring to Figs. 6 and, at step 601, initial values of 

grain, gmax, Gmin, G maX/ Pi and g ca i are set as follows. gmin is set 
as possible minimum gray level by selecting a lowest value of 
a gray level distribution curve on Fig. 4, and g max is set as 
possible maximum gray level by selecting a highest value of a 

25 gray level distribution curve on Fig. 4. And Gmax is set as 
equal to gmax and Gmin is set as equal to gmin- Also, Pi is set 
to mid point of g ma x and gmin as (gmax + gmin)/2. 



Furthermore, g ca i is set as equal to p ± . After 
initializing initial values, at step 602, entropy values 
E(gmin), E(g max ) and E(g ca i)of gmin, gmax and g ca i are computed. 
After computing entropy values E(g min ), E(g max ) and E(g ca i), 
5 E(g nii n) and E(g max ) are compared at step 603. 

If E (gmin) is higher than E(g ca i), the value of g t em P is 
changed to the value of g min and the value of g fix is set to G max 
at step 604. By changing values and not influencing to value 
of G max and Gmin/ new temporal values g t emp and g f i X are set as 
10 follows at step 604. 

If E(g min ) is equal or less than E(g ca i), the value of 
gtemp is changed to the value of g max and the value of g fix is 
set to G min at step 605. By changing values and not influencing 
to value of G max and Gmin, new temporal values g t em P and g fix are 
15 set as follows at step 605. 

After changing value of gtemp or g fix according to the 
comparison result and Pi is set to gmid/ gmid is computed by 
(gfix + gtemp) /2 and E m id is computed by (E(g te m P ) + E(g fix ))/2 at 
step 606. 

20 At step 607, p i+ i is computed by using a linear equation 

f with (gtemp, 0) and (gmid/ E (gmid) ) and E i+ i is set to E(p i+ i). The 
linear equation f is f(g)=ag+b. 

After computing Pi+i, it is compared with any two of 
previous Pi at step 608. 

25 At step 609, if there are identical two piS, it is ended, 

and at step 610, if there are not identical two PiS, g t em P is 
set to pi+i and g ca i is newly determined by (g t em P +gfix) /2, E(gmin) 
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is set to E i+ i and g t em P is set to P i+ i. After setting new value 
for g C ai, steps 602 and 608 are reputedly performed. For 
helping to understand steps for obtaining optimal threshold of 
Fig. 6, pseudo code is shown in below table. 



Table. 2 



Set flag=True; 






Set gmin= possible minimum occurring 


gray 


level; 


Set g m ax= possible maximum occurring 


gray 


level; 


Set G m i n =g m i n ; 






Set G m ax~gmaxJ 






Set pi=int [ (gmax+gmin) /2] ; 






set g C ai=Pi; 






Compute E (gmin) ; 






Compute E (g ma x) ; 






Compute E (g ca i) ; 






While (flage == True) 






if (E(g ca i)< E(g min )) then 






Set gtemp = gmiiw 






set gfix = G max / 






Else 






Set gtemp = gmax r 






set gf±x= Gmin; 






Set gmid = (gfix + gtemp) /2; 
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Set Pi=gmid; 

Set E^d = (E(g t emp)+E(g fix ) ) /2; 

Generate linear equation f using (gtemp, 0) and (gmidfEmid); 
Set p i+ i = f -1 (E(pi) ) ;Set E i+1 = E(p i+1 ); 

If (pi+i=any two previous p±) then 
set flag = false; 

Else 

set E (groin) = E i+i ; 

Set gtemp = Pi+lr gcal = ( gtemp + gfix)/2; 

END IF 
End While. 



As mentioned above, the present invention can quickly 
find the optimal threshold value by analyzing entropy 
characteristic of image based on a segmentation completion 
5 condition and a fixed point iteration. 

While the present invention has been described with 
respect to certain preferred embodiments, it will be apparent 
to those skilled in the art that various changes and 
modifications may be made without departing from the scope of 
10 the invention as defined in the following claims. 
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